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Abstract 



We use the scattering approach to investigate the nonhnear current-voltage 
characteristic of mesoscopic conductors. We discuss the leading nonlinearity 
by taking into account the self-consistent nonequilibrium potential. We em- 
phasize conservation of the overall charge and current which are connected to 
the invariance under a global voltage shift (gauge invariance). As examples, 
we discuss the rectification coefficient of a quantum point contact and the 
nonlinear current-voltage characteristic of a resonant level in a double barrier 
structure. 

PACS numbers: 72.20.Ht, 73.40.Ei, 73.40.Gk, 73.50.Fq 
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Typeset using REVT^ 



Introduction - The scattering approach is widely used to describe electric transport of 
phase-coherent mesoscopic conductors in the linear response regime Nonlinear effects 
of interest are asymmetric current- voltage characteristics and rectification [Q], the evolution 
of half-integer conductance plateaus 0, the breakdown of conductance quantization [Q, 
and negative differential conductance and hysteresis ||^. In general, in the nonlinear regime 
inelastic processes which destroy the phase coherence can play an important role. Below, 
however, we emphasize phase- coherent nonlinear transport. The nonequilibrium state is 
determined by the potential generated by the nonequilibrium charges piled up in a biased 
conductor. A reasonable theory of nonlinear electric transport has to take into account this 
potential self-consistently. Most importantly, for a mesoscopic sample there exists a Gauss 
volume which encloses the mesoscopic conductor and nearby gates (capacitors) such that the 
electric flux through the surface of this volume vanishes 0. Thus the total charge inside this 
volume is conserved. Consequently, the I-V-characteristic is gauge invariant: it is invariant 
under a global potential shift and depends only on the differences of the voltages applied 
to the contacts of the sample or nearby capacitors Nonlinearities in tunnel contacts 
have already been addressed in an early paper by Frenkel ||^ . The necessity to consider the 
self-consistent potential has been emphasized by Landauer [^]. Nevertheless, many works 
on nonlinear transport consider only the non-interacting case or consider interactions 
which are not gauge invariant. An exception is, e.g., a work by Kluksdahl et al. ^ which 
provides a self-consistent numeric treatment of a tunneling barrier with the help of Wigner 
functions. 

In this Letter we present a self-consistent, gauge-invariant theory of weakly nonlinear 
transport using the scattering approach. A gauge invariant theory of linear ac- conductance 
is the subject of Refs. P,p!OHr^. The discussion of non-linear transport presented here is 
based on similar concepts. We can illustrate the main points using an effective Hartree 
approach to treat the effects of the long range Coulomb interaction. The Hartree approach 
has a wide range of applicability but excludes single-charge effects [ITS . 



The leading nonlinearity - Consider a mesoscopic conductor which is connected via con- 



tacts a = 1, to electron reservoirs. We allow that some parts of the conductor are 
disconnected from other parts in order to include into the formalism the presence of nearby 
gates (capacitors). The transport properties are described by the scattering- matrix elements 
Sa(3nm. which relate the out-going current amplitude in channel n at contact a to the inci- 
dent current amplitude in channel m at contact f3. We denote by s^/s the scattering matrix 
with rows and lines associated with the channels n and m fl^. The scattering matrix is a 



function of the energy E of the wave and is a functional of the electric potential U{x, {V^7}) 
in the conductor. This potential in turn depends on the shifts eVy of the electrochemical 
potentials fi^ of the reservoirs away from the equilibrium state associated with /x^"^. Hence, 
the determination of the scattering matrix Sai3{E,{Vy}) as a function of the energy and 
the voltage shifts in the reservoirs defines a formidable self-consistent problem. Once this 
scattering matrix is found we can calculate the current through contact a 

2e ^ r 

Ia = -rY. dEf{E-EF- ^^pi^, {V^}) (1) 

where f{z) = (1 + exp{z / ksT))^^ is the Fermi function of a reservoir at temperature T, and 
AafsiE, {V^f}) = Tr[laSai3 — s^^(-E', {V^f})Saf3{E , {V.y})] are the screened (negative) transmis- 
sion functions. To discuss weakly nonlinear transport, we write for the current 

Jq, = ^ GapVp + ^ GafS-^fVjsV^ + ... . (2) 

The coefficients and Gap'y are obtained from an expansion of Eq. (|1|) with respect to 
the voltages V^. One obtains for the linear conductance the well-known expression Ga/s = 
{2e^ /h) J dE {—dEf)Aaj3, where the Aaji are taken at Vi = ... = Vn = 0. An expansion of 
Eq. dlD up to 0{Vf3V^) yields 

Ga^^ = jJdE i-dEf) {2dv,Aap + edEA^f35p^) . (3) 

Conservation of the total charge (gauge invariance) implies that the currents are conserved 



PJTm and that the currents be independent of a global voltage shift W: Va ^ Va + W . 
These conditions imply the sum rules |TT| Gap = H/3 Gap = 0. Similarly, the coefficients 



of the second-order nonlinearity must obey 




(4) 



The sum rules also follow from the unitarity of the scattering matrix, i.e. J2a = 
J2f3 ^a(3 = 0, and from the gauge invariance condition edEAap + J^-y dy^Aap = 0. 

To determine the derivatives dy^Aa/B we recall that the Aajj are functionals of the 
electric potential U{x,{V^}). Within linear response theory the potential variation 
AU (x) = U{x) — Ueq{x) away from equilibrium is determined by the characteristic potentials 
u^{x) = {dU{x)/dVy)eg. To derive the characteristic potentials we introduce the injectivities 
dn{x,a)/dE = —{l/Ani) J^isl^pai^^t^a/eSU^x)) — {6s^f^^/e6U{x))sfBoi]- The injectivity which is 
related to a dwell-time |]TO|,[T3| is a partial local density of states, i.e. the sum over all contacts 
a gives the local density of states dn{x)/dE. Variations eVa of the electrochemical potentials 
fia in the contacts at fixed electric potential Ueq{x) in the conductor induce a charge den- 
sity q^^\x) = e'^ Y.a{dn{x, a) / dE) V^. In general, however, the electric potential changes by 
AU{x) which induces an additional charge density q^'^\x) = —e'^JIl{x,y)AU{y)d^y, where 
n(x, y) is the Lindhard polarization function. Inserting the total nonequilibrium charge 
qW _|_ g(2) g^^d ^\^Q potential variation AU = UoYa into the Poisson equation yields for the 
characteristic potentials an equation with a nonlocal screening term 



For simphcity, the dielectric constant is to be unity. The derivative dy^Aap can be expressed 
in terms of the functional derivative 5Aaf3/5U{x) and the characteristic potentials, which 
yields for Eq. (H) 



Discrete potential model - We can obtain analytical results by using a convenient coarse 
graining model for the conductors We discretize the conductor by decomposing it 

into M regions where the electrostatic potential and the charge density are assumed to be 
constant. Space dependent quantities are written as M-dimensional vectors or as M x M- 
matrices. For example, the electrostatic potential variation AU{x) becomes a vector AU = 
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(AC/i, A[/m), where AUk denotes the electrostatic potential shift in region k. We write 
the injectivities as vectors = e^{dNia/dE, ...,dNMa/dE) where we multiplied by for 
later convenience. The functional derivative with respect to the potential U{x) is replaced 
by the gradient Vu, and the polarization function is a matrix 11. We introduce a geometric 
capacitance matrix C associated with these regions and which is calculated from the Poisson 
equation. The nonequilibrium charge distribution is q = CAU = J^p'DfsVfs — IIAU which 
yields for the characteristic potentials Ufs = (11 + C)~^D/3. Hence, Eq. (P) becomes 

2 

Ga(3y = X / dE{-dEf){2VuAap (H + C)-'B., - eSf^-ydEA^p) . (7) 
We mention that the charge qk in region k can be written in terms of an electrochemical 



capacitance matrix [T^, qk = Y^pC'^^Vp where = Y.iCkiUip. 

In the following we investigate two-terminal devices. Due to gauge invariance, the current 
= = GiiAF + Giii(A\/)2 + ... depends only on the voltage difference AV = Vi- V^. 
The second-order nonlinearity is thus determined by the single rectification coefficient Cm. 

Rectification of a quantum point contact - We follow closely Ref. and discretize the 



quantum point contact (QPC) by defining to the left and to the right of the constriction 
two regions Qi and ^2 with sizes of the order of the screening length. We assume zero 
temperature and take all quantities at the Fermi energy. A voltage difference across the 
QPC induces a dipole consisting of charges qi and q2 = —qi which reside in Qi and ^2, 
respectively. Each occupied subband j contributes with a transmission probability to 
the total transmission function T = Due to charge conservation, the geometrical 

capacitance matrix has the specific form Cm = (— 1)^^'C. Assuming a barrier potential of 
the constriction which allows a WKB approximation, we can neglect nonlocalities and write 
Uki = D^Ski- With each channel j one can associate a density of states Dl in region Q^, 
and it must hold Dk = J2j Dj. = J2a Dka- For a voltage variation Vi, channel j injects from 
contact 1 into Q2 a charge D^T^ /2)Vi and into Qi a charge -Di(l — /2)Vi. A sum over all 
channels j yields the injectivity Di of contact 1 with the components Du = Di{l — Ti/2) 
and D21 = D2T2/2. Here, Tk = D^' DiT^ 

are averaged transmission probabilities. We 
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obtain for the characteristic potentials uu = 1 — Ti/2 — C^/Di and U21 = T2/2 + j D2, 
where the electrochemical capacitance is given by = R/{C~^ + Df^ + D2^) and where 
we introduced the average reflection probability R = 1 — (Ti + T2)/2. The rectification 
coefficient becomes thus 

Cm = ^ ( cr(^^i^ - ^u.^) + [C"^(^2-' - D-,') - 2 - ^2)]e9^rj . (8) 

As one expects, Cm vanishes for a symmetric barrier. Asymmetries can be due to dif- 
ferences of Dki Tk, and duf,T for different k. Consider a single-channel conductor where 
Ti — T2 vanishes and which can be seen as a model for a one-dimensional conductor with 
an impurity. In the limit C ^ -D^, charges are screened and the rectification coefficient is 
Giu = {e^/h)R{du,T - du^T). On the other hand, in the limit Dk < C, Eq. (|) gives 
Cm = {e^/h){D2 — DiiROeT / [Di + D2) which is proportional to the asymmetry of the 
density of states. Clearly, with a discretization by more than two regions, an additional 
source of asymmetry can result from asymmetric electric screening (unequal capacitance 
coefficients). 

Resonant tunneling harrier - To model the resonant tunneling barrier, we consider three 
discrete potential regions denoted by indices 0, 1, and 2 which correspond to the well, the 
left side of the barriers, and the right side of the barriers, respectively. We assume fully 
screened charges. Consider a single resonant level with energy + cUq which is clearly 
separated from other resonant levels, from the band bottoms, and from the left and right 
barrier tops. The linear conductance of such a structure is given by Gu = {2e^/h)T where 
T is the transmission probability of the double barrier structure. Linear response is only 
valid if the bias |eAy| is much smaller than the width F of the resonance. To find the 
nonlinear characteristic, we restrict ourselves to near resonant conditions where the trans- 
mission probability can be approximated by the Breit-Wigner formula T{E, U) = rir2/|Ap 
with |Ap = {E — Er ~ gUq)^ + r^/4. In this range the nonlinearity is solely due to the 
resonance. In particular, we assume that Fi and and thus the width F = Fi + r2 of the 
resonant level and the asymmetry, AF = Fi — F2, are constant. The injectivities are |Tl| 



Doa{E,Uo) = (e^/27r)rQ,/| Ap. The relation between the electrostatic potential shift AUq 
and the voltage shifts = (/Xa — /i'^^)/e (= AUa) is obtained from the charge-neutrality 
condition 

D,,{E, Uo) dE + / Do2{E, Uo) dE - / D,{E, U^-^) dE = , (9) 

-oo J —oo J —oo 

which can be integrated analytically. Eq. (P) determines W = AUq — (Vi + V2)/2 which 
depends on AV only. Equation (|lD yields for the current I = Ii = —I2 

I = -^^^ (arctan[2(AE - e{W - AV/2))/T] - arctan[2(AE - e{W + AV/2))/T]) , 

(10) 

where AE = n^q — {gUq'^ + Er) is the equilibrium distance between the Fermi energy and the 
resonance. Without considering the self-consistent shift AUq one would get a wrong result 
which is not gauge invariant. The current given by Eq. ([lOD saturates at a maximum value 
proportional to 7r/2 — arctan(2Ai?/r). The conduction is optimal for AE = and Fi = F2 
when / = 2(e//;.)Farctan(eAy/F). In Fig. |I]we have plotted the characteristic as a function 
for an asymmetry AF/F = —1/3 and for various values of AE. For the case of complete 
screening considered here, the resonant level fioates up or down to keep the charge in the 
well constant. An expansion of the current yields Gm = —{e^/h){AT/T)dET which is in 
accordance with (^) and with Ref. |jl2| (thin dotted lines in Fig. |l]). The case of incomplete 
screening can similarly be treated with our approach. In this case, at large voltages, the 
resonance can eventually fall below the conductance band bottom of the injecting reservoir 
as is known from the semiconductor double barrier structures. In the general case, even an 
elastically symmetric resonance can be rectifying if the electrical screening is asymmetric. 

Conclusion - In this work we have emphasized that a reasonable theory of nonlinear 
electrical transport in mesoscopic conductors requires a self-consistent treatment of the 
long-range Coulomb interactions which insures gauge invariant results. We showed that 
the scattering approach to mesoscopic conduction provides an appropriate method if one 
introduces screened scattering matrices which depend on the self-consistent potential. 
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FIGURES 

FIG. 1. Current-voltage characteristic of a single level of the resonant tunneling barrier; 

Ti = 0.5meV, T2 = ImeV, AE = (solid), AE = -T (dashed); AE = -2r (dotted). Thin 
dotted lines indicate the quadratic approximation. 
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